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Thermal Model Additions 


Myers Water Film Model 
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Surface Water Shedding Model (calibrated) 
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Enhanced Evaporation 

- Chilton-Colburn analogy underestimates evaporation rate by 30% 
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Process for Comparison 


• Determine Internal Heat Transfer Coefficient from Dry 
Cases 


- All Cases Use Same Coefficients 

f \- 0.22 f \- 0.38 


Nu — 0.004 Re Pr 


UJ 
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• External Heat Transfer Coefficient is Forced Laminar 
Where There is No Ice 

• Run All Dry Cases To Ensure Correlation Matches 

• Run Wet Cases for Validation 
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Conditions Used For Thermal Comparison 


Case 

P(Pa) 

V(m/s) 

T(K) 

LWC 

MVD 

AOA 

t(min) 

Warm Hold(Ref) 

57295 

92.7 

264.5 

0.5 

20 

0 

7 

Warm Hold(Scale) 

98525 

54.3 

266.9 

0.85 

27.8 

0 

7 

Descent(Ref) 

69981 

92.7 

253.1 

0.15 

20 

0 

7 

Descent(Scale) 

97422 

66.8 

254.9 

0.21 

24.5 

0 

7 

Cold Hold(Ref) 

57295 

92.7 

247.4 

0.15 

20 

0 

10 

Cold Hold(Scale) 

98318 

54.6 

245.5 

0.25 

27.8 

0 

10 
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Warm Hold (Ref) - Dry 
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Warm Hold (ref) - Wet 


Warm Hold (ref) - wet 
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Warm Hold (Ref) Ice Shape Comparison 
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Warm Hold (Re Scale) - Dry 
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Warm Hold (Re Scale) - Wet 
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Warm Hold (Re Scale) Ice Shape 


• No Ice from Experiment nor from LEWICE 
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Temperature(K) 


Descent (Ref) 
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Descent (Ref) 
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Ice Growth in Impingement Direction 
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Temperature(K) 


Descent (Re Scale) 


Descent (scale) - wet 
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Descent (Re Scale) 
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Cold Hold (Ref) - wet 
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Cold Hold (Ref) 
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Temperature(K) 


Cold Hold (Re Scale) - wet 
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Cold Hold (Re Scale) 
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Observations from Thermal Analysis 


• Temperature Prediction is Very Good to Excellent for Most 
Cases 

• Warm Hold Cases Show Predicted Runback Ice Forward of 
Experiment 

- Peak Ice Thickness Higher for LEWICE 

• Descent and Cold Show Predicted Runback Ice Forms 
Slightly Behind Experiment 

- Peak Thickness Higher for LEWICE, Especially Upper Surface 

• Ice in Experiment Grows Toward Leading Edge While 
LEWICE always grows Ice Normal to Surface 

• Further Refinement of Runback Model May Be Necessary 

• External Heat Transfer Coefficients for Residual Ice Shapes 
Need to Be Separately Validated 
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Ice Breakup Model 


Breakup Threshold (Hauk) 
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Sticking Efficiency (Currie) 




m 


= 1 - 


^ COS (^))( 0 ' 57 + 7 ' 5 * l 0 ~ 4 [ V in,p C0S (^)l 5 


£ = - 0. 1425 + 47 . 292 TWC- 1 979 . 1 61 TWO 


- For TWC < 0.12 kg/m 3 and % = 0.14 for TWC > 0.12 
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Conditions for Ice Crystal Comparison 
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Ice Shape Prediction for Scan 996 
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Ice Shape Prediction for Scan 1003 
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Sticking Efficiency on Wedge at Various 

Particle Sizes 
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Ice Thickness (mm) 


Ice Thickness Prediction for Scan 796 

(NACA0012) 
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Observations from Ice Crystal Comparison 


• Peak Thickness is Over Predicted by LEWICE while Extent 
is Under Predicted 

- Additional Erosion Effects may be Needed 

- Improved Model for Reimpingement of Ice Crystals 

• Additional Data is Needed to Complete Model 
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